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1 Introduction 

A graph polynomial is an algebraic object associated with a graph that is usually 
invariant at least under graph isomorphism. As such, it encodes information 
about the graph, and enables algebraic methods for extracting this information. 
This chapter surveys a comprehensive, although not exhaustive, sampling of 
graph polynomials. It concludes Graph Polynomials and their Applications I: 
The Tutte Polynomial by continuing the goal of providing a brief overview of 
a variety of techniques defining a graph polynomial and then for decoding the 
combinatorial information it contains. 

The polynomials we discuss here are not generally specializations of the Tutte 
polynomial, but they are each in some way related to the Tutte polynomial, 
and often to one another. We emphasize these interrelations and explore how 
an understanding of one polynomial can guide research into others. Wc also 
discuss multivariablc generalizations of some of these polynomials and the theory 
facilitated by this. We conclude with two examples, one from biology and one 
from physics, that illustrate the applicability of graph polynomials in other 
fields. 

2 Formulating Graph Polynomials 

We have seen two methods for formulating a graph polynomial with the lin- 
ear recursion (deletion/contraction) and generating function definitions of the 
Tutte polynomial in the previous chapter. Here we will see several more. We 
begin with one of the earliest graph polynomials, the edge-difference polyno- 
mial, a multivariable polynomial defined as a product and originally studied 
by Sylvester [Syl78] and Peterson [Pet91] in the late 1800s. More recently, it 
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has been used to address list coloring questions (see Alon and Tarsi [AT92] and 
EUingham and Goddyn [EG96]), where a list coloring of a graph is a proper 
coloring of the vertices of a graph with the color of each vertex selected from a 
predetermined list of colors assigned to that vertex. 

Definition 1 The edge- difference polynomial. Let (ui,...,d„) he an ordering 
of the vertices of a graph G. Then D (G; xi, . . . , a;„) = Y\i<j i^i ~ ^j)' where 
the product is over all edges (vi,Vj) of G. 

Note that a proper coloring of G corresponds to finding positive integer 
values Ni (not necessarily distinct) for each of the Xi's so that D (G; A^i, . . . , Nn) 

There are also several polynomials based on various determinants (or even 
permanents; see Pathasarthy [Par89] for a survey) involving the adjacency ma- 
trix of a graph. Recall that A{G), the adjacency matrix of a graph, has entries 
Qij = 1 if is an edge of the graph and if it is not. The characteristic 
polynomial is the classic example of a such a graph polynomial, and will be 
discussed further in Section [S) 

Definition 2 The characteristic polynomial. Let A(G) he the adjacency matrix 
of a graph G. Then f{G\x)^\xI-A{G)\. 

Other examples of such polynomials are the idiosyncratic polynomial intro- 
duced by Tutte, see [Tut79], that is defined by u{G;x,y) = \A{G) + y{J - I - 
A) — xl\, where J is the matrix having all entries equal to 1. Also n{G; x, y) = 
\xI—D{G) + A{G)\ a polynomial introduced in [Kel65], where D{G) is the degree 
matrix of G that is the diagonal matrix with deg(?) in the position {i,i). Note 
that J — / — A is the adjacency matrix of the complement of G and when G is 
a simple graph, D{G) — A{G) is just the Laplacian matrix L{G) of the graph. 

A number of important graph polynomials may be defined by state model 
formulas. Loosely speaking, a state of a graph is some configuration resulting 
from making local assignments for substructures (e.g. the edges or vertices) of 
the graph. These assignments may be, for example, associating an element of a 
given set to each vertex, or even the result of reconfiguring the edges incident 
with a vertex. A graph polynomial is formed by associating an expression, often 
a weighted monomial, to each state of the graph, and then summing over all 
possible graph states. The language comes from physics, and is also found in 
knot theory. We will see several state model graph polynomials among those 
surveyed below, as well as an application of this method in the Potts model of 
statistical mechanics in Section [S] 

An early example of a graph polynomial given by a state model formulation 
is P(G; x), the Penrose polynomial. This polynomial graph invariant for planar 
graphs was defined implicitly by Penrose [Pen69] in the context of tensor dia- 
grams in physics, but an excellent graph theoretical exposition can be found in 
Aigner [Aig97]. To compute P{G;x), let G be a plane graph, and let Gm be 
its medial graph, face two-colored with the unbounded face colored white. At 
each vertex, we consider three possible local reconfigurations, as in Figure [l] 
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A state S of Gm then results from choosing one of these three reconfigurations 
at each vertex of Gm and consists of a set of disjoint closed curves (like a knot 
diagram). Furthermore, to each local reconfiguration at a vertex v, we assign a 
weight uj^S^v) that is +1, 0, or —1 for a white, black, or crossing configuration, 
respectively. 




Figure 1: The three possible local reconfigurations at a vertex v, identified, from 
left to right, as white, black, and crossing. Which strand passes over which in 
the crossing configuration does not affect the computation 



Definition 3 The Penrose polynomial. Let G he a planar graph with medial 
graph Gm, and let St{Gm) be the set of states of Gm and let St'{Gm) be the set 
of states with no black configurations. Then, 

St(G,„) \ \i>GG„ / / St'(G™) 

where k{S) is the number of components in the graph state S , and cr{S) is the 
number of crossing vertex configurations chosen in the state S . 

For example, if G is the 6'-graph consisting of two vertices joined by three 
edges in parallel, then P{G; x) — x^ — ix^ + 2x, as in Figure[2j The Penrose poly- 
nomial may also be computed via a linear recursion relation (see Jaeger [Jae90] 
for example). 

The Penrose polynomial has some surprising properties, particularly with 
respect to graph coloring. The Four Color Theorem is equivalent to showing 
that every planar, cubic, connected graph can be properly edge-colored with 
three colors. The Penrose polynomial, when applied to planar, cubic, connected 
graphs, encodes exactly this information (see Penrose [Pen69]): 



P(G;3) = 




P (G; —2) ~ the number of edge-3-colorings of G. 



3 Some Interrelated Polynomial Invariants 

We present a further sampling of graph polynomials here. They are each related 
in some way to the Tutte polynomial, and have additional relations among 
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Figure 2: Computing the Penrose polynomial of a graph G from the states of 
its medial graph 



themselves. These relations lead to combinatorial insights as results for any one 
polynomial then inform those related to it. 

3.1 Characteristic and Matching Polynomials 

The characteristic and matching polynomials are particularly interrelated, so we 
treat them together here, beginning with the characteristic polynomial /(G; x) 
already introduced in Definition [2] Note that /(G; x) is a monic polynomial of 
degree n. Furthermore, since the adjacency matrix A is real and symmetric, all 
its eigenvalues are real, and thus all the zeros of /(G; x) are real. 

By using properties of determinants we can find interpretations of the coef- 
ficients of /(G; x) in terms of the principal minors of A. A principal minor of 
order r is the determinant of an r x r submatrix of A obtained by choosing r 
rows and columns with the same set of indices. 

Proposition 1 Suppose that f{G;x) — X^iLo*^'^" Then (— l)*ai is equal to 
the sum of the principal minors of A with order i. 

This property of the characteristic polynomial can be found, for example, in 
Horn and Johnson [HJ90]. 

Since the diagonal elements of A are all zero, we have that ai = 0. The 
principal minors of order two and three which are not zero are of the form 
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I J — /|, where J is the matrix having all entries +1 and J — I has order 2 or 
3. The 2x2 submatrices J — / of A{G) correspond naturally to the edges of G 
and the 3x3 submatrices J — I correspond to the subgraphs of G. Thus 
C2 = —\E{G) \ and — C3 is twice the number of K^^ subgraphs of G. 

A linear subgraph of G is a subgraph whose components are edges or cycles. 
An expression for the coefficients of /(G; x) in terms of linear subgraphs is given 
in the following. 

Proposition 2 The coefficients of the characteristic polynomial may be ex- 
pressed as 

(-l)^a, = 5](-l)'-('^)2'-*(^), 

A 

where r is the rank function and the sum is over all linear subgraphs A of G 
having i vertices. 

Note that because A is a linear subgraph, r* (A) is simply the number of 
components in A that are cycles. The proof of Proposition |2] uses Proposition [T] 
and can be found in Harary [Har62], while a detailed history of this result is 
given by Cvetkovic, Doob, and Sachs [CDS80]. 

As with the Tutte polynomial we also have some reduction formulas and an 
expression for the derivative of the characteristic polynomial. 

Theorem 1 The characteristic polynomial of a graph satisfies the following 
identities: 

1. fiGUH;x) = fiG;x)f{H;x), 

2. f{G; x) — f{G \e,x) — f{G — u ~ v; x) if e = {u, v} is a cut-edge of G; 
3- fj{G;x)^j:.^y^^^nG-v;x). 

A proof of these properties can be found in Godsil [God93] . Ite m [T] is an 
easy exercises in matrix theory, as in Horn and Johnson [HJ90] . Item |2|can be 
proved by using Proposition [2] and considering the linear subgraphs of G that 
use the edge e and the ones that do not use it. The result follows because e is 
in no cycle of G if and only if it is a cut-edge. Item [3] can be proved by using 
Proposition [l] since any principal minor of order i is counted n — i times in the 
right hand side of the formula in Item [3j 

Given a graph G, the collection of (unlabeled) subgraphs G—v, for v e V{G), 
is called the deck of G, and the individual subgraphs are called cards. Thus, the 
deck for a graph on n vertices consists of n graphs, each of which has n — 1 ver- 
tices. Ulam's reconstruction conjecture in [Ula60] asserts that any finite graph 
G with more than two vertices is uniquely determined by its deck, see [Ula60], 
and we call any graph that satisfies the conjecture reconstructible. Similarly, an 
invariant of G which can be deduced from the deck is called reconstructible. 

Clearly, the number n of vertices of a graph is reconstructible. Also, as every 
edge is present in exactly n — 2 cards, the number of edges is also reconstructible. 
The following useful example of reconstructibility is due to Kelly [Kel57] . 
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Lemma 1 (Kelly's lemma) Let G and H be graphs, and let v{H, G) denote 
the number of subgraphs of G isomorphic to H. Then 

{\ViG)\-\ViH)\HH,G)= J2 HH,G-v). 

v<£V{G) 

The proof is by a double counting argument, and it follows that v{H,G) is 
reconstructible whenever < |y(G)|. 

Tutte proved in [Tut67] that the Tutte polynomial is reconstructible, and 
thus the chromatic polynomial, the flow polynomial, the number of spanning 
trees and any invariant mentioned in the previous chapter are also reconstructible. 
Tutte also proved that the characteristic polynomial is reconstructible in [Tut79] . 

Theorem 2 The characteristic polynomial of a graph is reconstructible. 

For the proof, note that we have immediately from Theorem [l] that f'{G;x) 
is reconstructible. It then remains to prove that the constant term of /(G; x) is 
reconstructible. But by Proposition [T| this is the same as proving that |A(G)| 
is reconstructible. Then, using Theorem [2] and an extension of Kelly's lemma 
(see [Koc81]), the problem is reduced to proving that the number of Hamiltonian 
cycles is reconstructible. A complete proof of Theorem [2] based on the proof in 
Kocay [Koc81], can be found in [God93]. 

Let us turn now to the matching polynomial. An i-matching in a graph G 
is a set of i edges, no two of which have a vertex in common. Let ^i{G) denote 
the number of i-matchings, and set ^q(G) = 1. Thus <i>i(G) = m is the number 
of edges of G, and if n, the number of vertices, is even, then <i>„/2(G) is the 
number of perfect matchings of G. 

Definition 4 Let G be a graph. Then the matching polynomial of G is 

^iiG■,x) = J2i-^yMG)x--'\ 

i>0 

A more natural polynomial might be the matching generating polynomial, 
given as the generating function of i-matchings by 

g(G;a;)=^<i>,(G)x\ 

i>0 

However, the two polynomials are related by the identity 

/.(G;x) = a;"5(G; {-xf), 

so there is no essential difference between them. 

The matching polynomial is also known as the acyclic polynomial in Gut- 
man and Trinajstic [GT76], matching defect polynomial in Lovasz and Plum- 
mer [LP86] and reference polynomial in Aihara [Aih76]. It has appeared in- 
dependently in several different contexts. In combinatorics, it was probably 
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introduced by Farrell in [Far79a], but since the matching polynomial is essen- 
tially the same as the rook polynomial for bipartite graphs (see Farrell [Far88]), 
then its origin can be traced back at least to Riordan [Rio58]. In statistical 
physics it appears because of the monomer-dimer problem and was introduced 
by Heilmann and Lieb in [HL70] and independently by Kunz in [Kun70] . Finally, 
in theoretical chemistry was introduced by Hosaya in [Hos71] and later in con- 
nection with the so-called topological resonance energy by Gutman, Milun and 
Trinajstic in [GT76, GMT76, GMT77] and independently by Aihara in [Aih76]. 
For a full account of the history of the matching polynomial see Gutman [Gut91]. 

As with the Tutte and characteristic polynomials, we have some reduction 
formulas for the matching polynomial. The proof of the following theorem can 
be found in Godsil [God93]. 

Theorem 3 The matching polynomial satisfies the following identities: 

1. fi{GU H;x) = ^i{G;x)^{H;x), 

2. fJ.{G; x) = /i(G \e;x) ~ ^{G — u — v;x) if e = {u, v} is an edge of G, 

3. fi{G; x) = xfj.{G -u;x) - J2{u,v}eE{G) -v-u;x), if u e V{G), 

4. ^tJ-iG;x)=J2vev{G)l^(C!~v;x). 

To choose an i-matching in GU H you need to choose an s-matching in G and 
a i-matching in H such that s + t = i. Item [T] then follows by the fundamental 
counting principle. For Item [2] notice that the set of i-matchings can be parti- 
tioned into those i-matchings that use the edge e and those that do not use it. 
Item [3] follows similarly. Finally, every z-matching of G with z < n/2 is counted 
n — 2i times in the right-hand side of the formula in Itcm|4] so the result follows. 

When the graph G is a forest, a linear subgraph of G with j vertices cor- 
responds to a matching covering j vertices, with j even. Thus, Proposition [2] 
has the following corollary, observed by Hosaya [Hos71] and by Heilmann and 
Lieb [HL72]. 

Corollary 1 If G is a forest then f{G; x) = IJ.{G; x). 

An unexpected property of the matching polynomial, proved by Heilmann 
and Lieb [HL72], is that all its zeros are real, and furthermore the zeros for any 
graph G interlace with the zeros of any of the cards in its deck. The same paper 
also gives bounds for the zeros. 

Theorem 4 For any graph G, the matching polynomial /i(G; x) has only real 
zeros. Furthermore, if u is any vertex in G and i/ai, 02, . . . , a„ are the zeros of 
fJ-{G; x) while the zeros of ^{G — u; x) are a[, a'2, ■ ■ ■ , a'^-i; then 

0-1 l£ l£ 0,2 l£ 0,2 — ■ ■ ■ — '^n-l — '^ni 

that is, the zeros of ijl{G]x) and fi{G — u;x) interlace. 
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Theorem 5 The (real) zeros, ai, 02, . . . , a„, of ^{G;x), satisfy 

\ai\ < 2yJ maxdeg{G) — 1. 

We outline a proof of Theorem |4] from Godsil [GodSlb] that uses some of 
the results already mentioned for /i(G; x) and f(G;x). First, given a graph G 
and a vertex u in G, the path tree T{G, u) is the tree that has as its vertices 
the paths in G which start at u, and where two such vertices are joined by an 
edge if one represents a maximal proper subpath (i.e. all but the last edge) of 
the other. We then have the following proposition from [GodSlb] that leads to 
a proof of Theorem |4] 

Proposition 3 Let u he a vertex in a graph G, let T — T(G, u) be the path tree 
of G with respect to u, and let u' he the vertex of T corresponding to the path of 
length beginning at u. Then 

li{G — u]x) ^{T—u';x) f{T — u'\x) 
KG;x) " fi{T;x) " f{T;x) 

The last equality follows from Corollary [l] Because all the roots of the charac- 
teristic polynomial are real, we conclude that all zeros and poles of the rational 
function fiiG — u; x) / ^{G\ x) are real. An induction argument on the number 
of vertices in G then yields the conclusion that all the zeros of ii{G;x) are real. 

An interesting combinatorial consequence of Theorem |4] is the following re- 
sult of Heilmann and Lieb [HL72] , which gives a stark contrast with how little 
is known about the coefficients of the chromatic polynomial. 

Theorem 6 For any graph G, the sequence <i>o(G'), <i>i(G'), ■■■of coefficients 
of g{G;x) is log-concave, that is <i>f > <i>i_i<i>j_(-i. 

The characteristic polynomial has been well studied, particularly with re- 
spect to graphs with the same characteristic polynomial. Godsil [God93] gives 
a thorough treatment of both the characteristic and the matching polynomials. 
Another good reference for the characteristic polynomial is Biggs [Big96]. Just 
as the matching polynomial is a way to study the matchings of a graph, the 
characteristic polynomial is a way to study the spectra of the adjacency matrix 
of a graph. Cvetkovic, Doob, and Sachs have written a book [CDS80] dedi- 
cated to the spectra of the adjacency matrix, and Lovasz and Plummer [LP86] 
have a book devoted to the theory of matchings. Furthermore, although the 
characteristic polynomial is not a complete invariant of graphs, it is conjectured 
that the characteristic polynomial of a graph G is reconstructible from its poly- 
nomial deck, i.e. from the set of characteristic polynomials of the cards of G. 
See Gutman and Cvetkovic [GC75] for the conjecture, and then Cvetkovic and 
Lepovic [CL98], where it is proved in the case of trees. 

3.2 Ehrhart Polynomial 

A convex polytope P is the convex hull of a finite set of points in M™. We denote 
the interior of P (in the usual topological sense) by P°. A convex polytope P 
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is said to be a rational or integral polytope if all its vertices have rational or 
integral coordinates, respectively. We write d = dim P and call P a d-polytope. 

For P c a rational rf-polytope and t a nonnegative integer we define the 
functions i{P\ t) = \tP n Z™| and i{P; t) = \tP° n Z™|, where tP ^ {ta\a G P] 
is the i-fold dilatation of P. Ehrhart proved in [Ehr67a, Ehr67b] that these 
functions are quasi-polynomials, that is, they are of the form 

Cd{t)t'^ + ct-i{t)t'^-^ + ... + co{t), 

where each Ci{t) is a periodic function with integer period. Since i{P;t) is a 
quasi-polynomial, it can be defined for all t € Z. In fact, we have the following 
reciprocity law due to Ehrhart [Ehr67c]: 

i{P;-t) = {-l)n{P;t). 

For more on this beautiful theory see the monograph by Ehrhart [Ehr77]. 

From Ehrhart [Ehr67a,Ehr67b] (also see Stanley [Sta96]), we have that when 
P C is an integral c?-polytope, then i{P; t) and i{P; t) are polynomials, which 
leads to the following definition of the Ehrhart polynomial. 

Definition 5 Let P he an integral convex d-polytope. Then the Ehrhart poly- 
nomial of P is 

l{P] t)=Cn+Clt+ ... + Crf_lt"-1 + Cdt". 

From the early works of Ehrhart [Ehr67a] and Macdonald [Mac71] it is known 
that Co = 1 and Cd = vol(P), and that Cd-i is half of the surface area of 
P, normalized with respect to the sub-lattice on each face of P. Specifically, 
Cd-i = 1/2 void- where F ranges over all facets of P and the volume 
of a facet is measured intrinsically with respect to the lattice Z"* n Lp, where 
Li? is the afhne hull of F. The other coefficients were not well understood, until 
the later work of Betke and Kneser [BK85], Pommersheim [Pom93], Kantor 
and Khovanskii [KK93] and Diaz and Robins [DR97], but such interpretations 
go beyond the scope of this chapter. For the complexity of computing these 
coefficients see Barvinok [Bar94]. 

In the special case that P is a zonotope there is a combinatorial interpretation 
for the coefficients of the Ehrhart polynomial. First recall that if ^ is an r x to 
real matrix written in the form A= [ai, a^.], then it defines a zonotope Z{A) 
which consists of those points p of R'" which can be expressed in the form 

m 

p = ^Aiaj, < Ai < 1. 

1=1 

In other words, Z(A) is the Minkowski sum of the line segments [0,0^], for 
1 < i < n. For more on zonotopes, see McMuUen [McM71]. 

When A has integer entries, Stanley [StaSO], using techniques from Shep- 
hard [She74], proved that i{P;t) = Ex/(^)^'^'' where X ranges over all lin- 
early independent subsets of columns of A and where f{X) denotes the greatest 
common divisor of all minors of sizes |^| of the matrix A. 
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When A is a totally unimodular matrix, that is, the determinant of every 
square submatrix is or ±1, then Z{A) is described as a unimodular zonotope. 
For these polytopes the previous result shows that 

r 

fc=0 

where fk is the number of subsets of columns of the matrix A which are linearly 
independent and have cardinality fc. In other words, the Ehrhart polynomial 
i{Z(A);t) is the generating function of the number of independent sets in the 
regular matroid M{A). 

The incidence matrix DiG) of a graph G is totally unimodular, a long- 
standing result due to Poincare [PoiOl] with a modern treatment given by 
Biggs [Big96] . A linearly independent subset of columns in D corresponds to a 
subset of edges with no cycle. Thus, the coefhcient fk in this case is the number 
of spanning forests of G with exactly k edges. From the previous chapter we 
know that T{G] x + l,y) — X]fc=o fkX^~''t where r is the rank of the graph G. 
With these ingredients we get the following relations with the Tutte polynomial 
from Welsh [Wel97]. 

Theorem 7 If G is a graph and D is its incident matrix then the Ehrhart 
polynomial of the unimodular zonotope Z{D) is given by 

*(Z(Z?);t)=rT(G; 1+^,1), 

where r is the rank of G. 

In this case, the zonotope Z{D) is a r-polytope in R", where n is the number 
of vertices of G. 

The reciprocity law of Theorem [7] leads to the following geometric result, 
also from Welsh [Wel97]. 

Corollary 2 If D is the incidence matrix of a rank r graph G with n vertices 
then for any positive integer A the number of lattice points of lying strictly 
inside the zonotope tZ{D) is given by 

l{Z{D);t) = {-tYT{G-l-\,l). 

In particular we have that the number of lattice points strictly inside Z(D) is 
(^1)'^T(G;0,1). 

3.3 The Topological Tutte Polynomial of Bollobas and Ri- 
ordan 

The classical Tutte polynomial discussed in the previous chapter is an invariant 
of abstract graphs, so it encodes no information specific to graphs embedded 
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in surfaces. In [BR01,BR02], BoUobas and Riordan generalize the classical 
Tutte polynomial to topological graphs, that is, graphs embedded in surfaces. 
In [BROl], BoUobas and Riordan define the cyclic graph polynomial, a three 
variable deletion/contraction invariant for graphs embedded in oriented sur- 
faces. They extend this work in [BR02], using a different approach, with the 
four variable ribbon graph polynomial. Both of these polynomials extend the 
classical Tutte polynomial, but in such a way that topological information about 
the embedding is encoded. The version for oriented surfaces is subsumed by the 
version for arbitrary surfaces, so we focus on the latter here. The ribbon graph 
polynomial is also sometimes called the BoUobas- Riordan polynomial after the 
authors or the topological Tutte polynomial to emphasize that it simultaneously 
encodes topological information while generalizing the classical Tutte polyno- 
mial. 

First recall that a cellular embedding of a graph in an orientable or unori- 
entable surface can be specified by providing a sign for each edge and a rotation 
scheme for the set of half edges at each vertex, where a rotation scheme is simply 
a cyclic ordering of the half edges about a vertex. This is equivalent to a ribbon 
(or fat) graph, which is a surface with boundary where the vertices are repre- 
sented by a set of disks and the edges by ribbons, with the ribbon of an edge 
with a negative sign having a half- twist. This can also be thought of as taking 
a shght 'fattening' of the edges of the graph as it is embedded in the surface, or 
equivalently as 'cutting out' the graph together with a small neighborhood of 
it from the surface. Figure |3] shows a graph with two vertices and two parallel 
edges, one positive and one negative. It is embedded on a Klein bottle, and the 
ribbon graph is a Mobius band with boundary. 



Figure 3: A ribbon graph which is a Mobius band with boundary 

In addition to the usual graphic characteristics such as number of vertices, 
connected components, rank, and nullity for a ribbon graph G wc also consider 
be (G), the number of boundary components of the surface, and t(G), an index of 
the orientability of the surface. The value of t{G) is if the surface is orientable, 
and 1 if it is not. Thus, t{G) is 1 if and only if for some cycle in G, the product 
of the signs of the edges is negative. 

Definition 6 Let G be a ribbon graph, that is, a graph embedded in a surface. 




11 



Then the topological Tutte polynomial of Bollobds and Riordan is given by 
R{G; X, y, z,w) = Yl " 

AQE{G) 

e Z[x,y,z,w]/{w'^ - w). 

As previously, r{A), k{A), n{A), and now also bc{A) and t{A), refer to the 
spanning subgraph of G with edge set A, here with its embedding inherited from 
G. 

Clearly, by comparing with the rank-nullity generating function definition of 
the classical Tutte polynomial given in the previous chapter, this generalizes the 
classical Tutte polynomial. Like the classical Tutte polynomial, R{G; x, y, z, w) 
is multiplicative on disjoint unions and one point joins of ribbon graphs. More 
importantly, it retains the essential properly of obeying a deletion/contraction 
reduction relation. 

We must first define deletion and contraction in the context of embedded 
graphs. The ribbon graph resulting from deleting an edge is clear, but contrac- 
tion requires some care. Let e be a non-loop edge. First assume the sign of e 
is positive, by flipping one cndpoint if necessary to remove the half twist (this 
reverses the cyclic order of the half edges at that vertex and toggles their signs). 
Then G/e is formed by deleting e and identifying its endpoints into a single 
vertex V. The cyclic order of edges at v comes from the original cyclic order at 
one endpoint, beginning where e had been, and continuing with the cyclic order 
at the other endpoint, again beginning where e had been. 

Theorem 8 If G is a ribbon graph, then 

R{G; X, y, z, w) = R{G/ e; x,y,z,w) + R{G - e;x,y,z,w) 

ife is an ordinary edge and R{G; x, y, z, w) = xR{G/e; x, y, z, w) if e is a bridge. 

The proof depends on a careful analysis of how each of the relevant parame- 
ters r{A), k{A), n{A), bc{A) and t{A) changes with the deletion or contraction 
of an edge. 

Repeated application of this theorem reduces a ribbon graph to a disjoint 
union of embedded blossom graphs, that is, graphs each consisting of a sin- 
gle vertex with some number of loops. Because of the embedding, the loops 
are signed, and there is a rotation system of half-edges about the single ver- 
tex. Not surprisingly, the topological information is distilled into these minors 
of the original graph, and to complete a deletion/contraction linear recursion 
computation, it is necessary to specify an evaluation of these terminal forms. 

Signed chord diagrams provide a useful device for determining the relevant 
parameters of an embedded blossom graph. Recall that a chord diagram consists 
of a circle with n symbols on its perimeter, with each symbol appearing twice 
and a chord drawn between each pair of like symbols. A signed chord diagram 
simply has a sign on each chord. A signed chord diagram D corresponds to an 
embedded blossom graph G by assigning a symbol to each loop and arranging 
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them on the perimeter of the circle in the chord diagram in the same order as 
the cychc order of the half-edges about the vertex. A chord receives the same 
sign as the loop it represents. If we 'fatten' the chords as in Figure |4] with 
a negative chord receiving a half-twist, then bc{G) is equal to the number of 
components in the resulting diagram, which is denoted hc{D). Similarly, since 
G has only one vertex, n{G) is the number of edges of G, which is the number 
of chords of D, so we denote this by n{D). We also set t{D) = t{G), and note 
that t{D) = if all chords of D have a positive sign, and t{D) = 1 otherwise. 
This, combined with the definition of R{G; a;, y, z, w) above, gives the following 
evaluation for these terminal forms. 




Figure 4: A signed blossom graph, its signed chord diagram, and the boundary 
components of the signed chord diagram 



Theorem 9 IfGis an embedded blossom graph with corresponding signed chord 
diagram D, then 

R{G;X,y,Z,w)^ ^ y^iD') ^J-bciD')+niD')^t(D') ^ 
D'CD 

where the sum is over all subdiagrams D' of D. 

Theorems [s] and [9] taken together give a linear recursion definition for R{G; 
X, y, z, w). There are a number of technical considerations, similar to the care 
that must be taken in contracting edges, but nevertheless many other prop- 
erties analogous to those of the classical Tutte polynomial hold. For exam- 
ple, R{G; X, y, z, w) has a spanning tree expansion, a universality property, and 
duality relation (in addition to the Bollobas and Riordan's originating work 
in [BR01,BR02], see also recent work by Chmutov [Chm] and Moffatt [Mof08]). 
Furthermore, Chmutov and Pak [CP07], and Moffatt [Mof, Mof08] have shown 
that R{G;x,y,z,w) also extends the relation between the classical Tutte poly- 
nomial and the Kauffman bracket and Jones polynomial of knot theory due to 
Thistlcthwaitc [Thi87] and Kauffman [Kau89]. 

3.4 Martin, or Circuit Partition, Polynomials 

In his 1977 thesis, Martin [Mar77] recursively defined polynomials M{G, x) 
and m{G; x) that encode, respectively, information about the families of cir- 
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cuits in 4-rcgular Eulcrian graphs and digraphs. Las Vcrgnas subsequently 
found a state model expression for these polynomials, extended their properties 
to general Eulerian graphs and digraphs, and further developed their theory 
(sec [Las79,Las83,Las88]). Both Martin [Mar 78] and Las Vcrgnas [Las88] found 
combinatorial interpretations for some small integer evaluations of the polyno- 
mials, while combinatorial interpretations for all integer values as well as some 
derivatives were given in [E-M00,E-M04a,E-M04b], and by BoIIobas [BoI02]. 

Transforms of the Martin polynomials, J{G;x) and j{G;x), given in [E- 
M98], and then aptly named circuit partition polynomials in [ABSOO], facilitate 
these computations, and for this reason wc give the definitions below in terms 
of J and j. Like many of the polynomials surveyed here, the circuit partition 
polynomials have several definitions, including linear recursion formulations, 
generating function formulations, and state model formulations. Wc give the 
state model definition, and refer the reader to [E-M04a, E-M04b] for the others. 

As with other state model formulations, we must first specify what we mean 
by a state of a graph (or digraph) in this context. Here an Eulcrian graph must 
have vertices all of even degree, but it need not be connected. An Eulerian 
digraph must have the indegree equal to the outdegree at each vertex, and 
again need not be connected. 

Definition 7 An Eulerian graph state of a/n Eulerian graph G is the result of 
replaeing each 2n-valent vertex v of G with n 2-valent vertices joining pairs 
of half edges originally adjacent to v. An Eulerian graph state of an Eulerian 
digraph G is defined similarly, except here each incoming half edge must be paired 
with an outgoing edge. 

Note that a Eulerian graph state is a disjoint union of cycles, each consis- 
tently oriented in the case of a digraph. 

Definition 8 The circuit partition polynomial. Let G be an Eulerian graph, 
let St{G) be the set of states of G, and let c{S) be the number of components 
in a state S G St{G). Then the circuit partition polynomial has a state model 
formulation given by 

J{G;x)= 

seSt{G) 

The circuit partition polynomial is defined similarly for Eulerian digraphs as 
j{G;x)= 

sest{G) 

The transforms between the circuit partition polynomials and the original 
Martin polynomial, as extended to general Eulerian graphs and digraphs by Las 
Vergnas, are: 

J{G; x) = xM (G; a: + 2) , for G an Eulcrian graph, and (3.1) 
j(G; x) = xm (G; a; + 1) for G an Eulerian digraph. (3.2) 
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The circuit partition polynomials have 'splitting' formulas, analogous to 
Tutte's identity for the chromatic polynomial given in the previous chapter, 
proofs for which may be found in [E-M98,E-M04b]. These formulas derive from 
the Hopf algebra structures of the generalized transition polynomial discussed in 



Subsection 4.3 but may also be proved combinatorially, as by BoUobas [Bol02] 



and [E-M04b]. 

Theorem 10 Let G be an Eulerian graph and G he an Eulerian digraph. Then 

J(G;x + j/) = ^J(^;a;) J(A^;y), 

where the sum is over all subsets A C E{G) such that G restricted to both A 
and A'^ = E{G) ^ A is Eulerian. Also, 

j{G; x + y) = J2j j (^'; y) > 

where the sum is over all subsets A C E{G) such that G restricted to both A 
and A'^ is an Eulerian digraph. 

The connection between the circuit partition polynomial of a digraph and 
the Tutte polynomial of a planar graph G is through the oriented medial graph 
Gm described in the previous chapter. Martin [Mar 77] proved the following. 



which we extend to the circuit partition polynomial via (3.2 1 



Theorem 11 Let G be a connected planar graph, and let Gm be its oriented 
medial graph. Then relationships among the Martin polynomial, circuit partition 
polynomial, and Tutte polynomial are: 

j{G; x) = xm{Gm', a; + 1) = xt{G\ x + l,x + 1). 

The proof of this theorem depends on a fundamental observation relating 
deletion/contraction in G with choices of configurations at a ve rtex in an Eule- 
rian graph state of Gm, as illustrated in Figure [s] Theorems 10 and [ll] combine 



to give the basis for many of the combinatorial interpretation of the Tutte poly- 
nomial along the line y = x described in the previous chapter. For more details, 
see Martin [Mar77,Mar78], Las Vergnas [Las79,Las83,Las88], BoUobas [Bol02], 
and also [E-M98,E-M00,E-M04a,E-M04b]. 

Evolving from the relation between the Tutte and Martin polynomials is the 
theory of isotropic systems, which unifies essential properties of 4-regular graphs 
and pairs of dual binary matroids. A series of papers throughout the 1980's 
and 1990's, including work by Bouchet [Bou87a], [Bou87b], [Bou88], [Bou89], 
[Bou91], [Bou93], as well as Bouchet and Ghier [BG96], and Jackson [Jac91], 
significantly extends the relationship between the Tutte polynomial of a planar 
graph and the Martin polynomial of its medial graph via the theory of isotropic 
systems. 
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Figure 5: An edge e in a planar graph G, with the corresponding vertex v in 
the oriented medial graph Gm (dotted edges). Deleting e corresponds to one 
possible configuration at v in an Eulerian graph state of Gm, while contracting 
e corresponds to the other 



3.5 Interlace Polynomial 

In [ABSOO], Arratia, BoUobas and Sorkin defined a one- variable graph poly- 
nomial motivated by questions arising from DNA sequencing by hybridization 
addressed by Arratia, BoUobas, Coppersmith and Sorkin in [ABCSOO], an appli- 
cation we will return to in Sectionjsj In [ABS04b], Arratia, BoUobas, and Sorkin 
defined a two-variable interlace polynomial, and showed that the original poly- 
nomial of [ABSOO] is a specialization of it, renaming the original one-variable 
polynomial as the vertex-nullity interlace polynomial due to its relationship 
with the two-variable generalization. 

Remarkably, despite very different terminologies, motivations and approaches, 
the original vertex-nullity interlace polynomial of a graph may be realized as 
the Tutte-Martin polynomial of an associated isotropic system (see Bouchet 
[Bou05]). For exploration of this relationship, see the works mentioned in Sub- 
section |3.4[ as well as Aigner [AigOO] , Aigner and Mielke [AMOO] , Aigner and 
van der Hoist [AvdII04], AUys [A1194], and also Bouchet 's series on multima- 
troids [Bou97, Bou98a, Bou98b, BouOl] . 

Both the vertex-nullity interlace polynomial of a graph and the two- variable 
interlace polynomial may be defined recursively via a pivot operation. This 
pivot is defined as follows. Let vw be an edge of a graph G, and let A^, 
and Ayuj be the sets of vertices in T^(G) \ {Vjw} adjacent to v only, w only, 
and to both v and w, respectively. The pivot operation "toggles" the edges 
among A^ , A^ and A^^ , by deleting existing edges and inserting edges between 
previously non-adjacent vertices. The result of this operation is denoted G""". 
More formally, G™ has the same vertex set as G, and edge set equal to the 
symmetric difference E{G)AS, where S is the complete tripartite graph with 
vertex classes A^ , A^ and Ayiu . See Figure [6] 
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Figure 6: Pivoting on the edge vw. Ay, A^, and A.^.^, are the sets of vertices 
of G adjacent to v only, w only, and to both v and w, respectively. These sets 
are constant in all the diagrams. Vertices of G adjacent to neither v nor w 
are omitted. Heavy lines indicate that all edges are present, and dotted lines 
represent non-edges. Note interchange of edges and non-edges among Ay, Ayj 
and A-ynj 



Also, G" is the local complementation of G, defined as follows. Let N(a) be 
the neighbors of a, that is, the set {w £ V : a and w are joined by an edge}. 
The graph G" is equal to G except that we "toggle" the edges among the 
neighbors of a, switching edges to non-edges and vice-versa. 

Definition 9 Let G be a graph of order n, which may have loops, hut no multi- 
ple loops or m,ultiple edges. The two-variable interlace polynomial may be given 
recursively by q{En) = y" for En, the edgeless graph on n > vertices, with 

q(G) = q{G -a)+ qiG'^' - b) + ((x - l)^ - l)q{G''>' -a-h), 

for any edge ah where neither a nor h has a loop, and 

q{G) = q{G-a) + {x-l)q{G--a), 

for any looped vertex a. 

Alternatively, the interlace polynomial has the following generating function 
representation. 

Definition 10 Let G be a graph of order n, which may have loops, but no 
multiple loops or multiple edges. Then the two-variable interlace polynomial 
may be given by 

q{G;x,y)= ^ (a; - If^^l-Hy - 
Scy(G) 

where r{G\g) and n{G\g) = \S\ — r{G\g) are, respectively, the ¥2-rank and 
nullity of the adjacency matrix of G\g, the subgraph of G restricted to S. 
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Definition 11 The vertex-nullity interlace polynomial is defined recursively as: 

\ x" if G = En, the edgeless graph on n vertices 
\qN{G~v;x) + qNiG'''^ -w;x) ifvweE{G). 

This polynomial was shown to be well defined by Arratia, BoUobas, and 
Sorkin for all simple graphs in [ABSOO], and then was shown in [ABS04b] to be 
a specialization of the two-variable interlace polynomial as follows. 

qN{G;y)^q{G;2,y)^ ^ (y-ir(«l»-). 

WCV{G) 

An equivalent formulation for qNiG; x) for simple graphs is given by Aigner 
and van der Hoist in [AvdH04]. 

A somewhat circuitous route through the circuit partition polynomial relates 
the vertex-nullity interlace polynomial to the Tutte polynomial. First recall 
that a circle graph on n vertices is a graph G derived from a chord diagram. 
Two vertices v and w in G share an edge if and only if their corresponding 
chords intersect in the chord diagram. Note that G is necessarily simple. 

For circle graphs, the vertex- nullity interlace polynomial and the circuit par- 
tition polynomial are related by the following theorem, noting that although G 
may be a multigraph, H is necessarily simple. 

Theorem 12 (Arratia, Bollobds and Sorkin [ABSOO], Theorem 6.1). If G is a 
A-regular Eulerian digraph, C is any Eulerian circuit of G, and H is the circle 
graph of the chord diagram determined by C , then j{G; x) = xq]^(H; x + I). 

This now allows us to relate the vertex-nullity interlace polynomial to the 
Tutte polynomial, a relation proved in [E-MS07] and also observed by Arratia, 
BoUobas and Sorkin at the end of Section 7 in [ABS04b]. 

Theorem 13 If G is a planar graph, and H is the circle graph of some Eulerian 
circuit of Gm, then q]^{H;x) — t{G;x,x). 

Proof. By Theorem [l2j j{Gm',x) = xqN{II; x + 1), but recalling that the cir- 
cuit partition and Mart in polynomials are simple translations of each other, we 
have from Theorem [Ulth^t /(G™; X) = xm(G„; x + 1), and hence ; .) = 

m{Gm',x) = t{G;x,x). 

The interlace polynomial has generated further interest and other applica- 
tions in Balister, BoUobas, Cutler, and Pebody [BBCP02], and Balister, Bol- 
lobas, Riordan, and Scott [BBRSOl], Glantz and Pehllo [GP06], EUis-Monaghan 
and Sarmiento [E-MS07]. 

4 Multivariable Extensions 

Multivariable extensions have proved valuable theoretical tools for many of the 
polynomials we have seen since they capture information not encoded by the 
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original polynomial. More critically, powerful algebraic tools not applicable to 
the original polynomial may be available to the multivariable version, provid- 
ing new means of extracting combinatorial information from the polynomial. 
While the multivariable indexing may make the defining notation somewhat 
bulky, these generalizations are natural extensions of classical versions, com- 
puted in exactly the same ways, only now also keeping track of some additional 
parameters in the computation processes. 



4.1 Generalized coloring polynomials 
and the U-polynomial 

The evaluation of the chromatic polynomial at A can be written as 

Xg(A)= 1- (4.1) 

0:V^{1,...,A} 
proper 

This was generalized to a symmetric function over (commuting) indetermi- 
nates Xi,X2, - ■ ■ by Stanley [Sta95] in the following way. 

Definition 12 Let G = {V, G) he a graph, let cf) : V V = {1,2, . . .}, and de- 
note the product Jl^ey ^<P{v) Then the symmetric function generalization 
of the chromatic polynomial is 



proper 



That this is a generalization of the chromatic polynomial can be seen by 
setting Xi — I for 1 < i < A and Xj = for j > A and noting that the expression 



in (4.1 1 for the chromatic polynomial evaluated at A results. 

Generalizing polynomial graph invariants is not a theoretical exercise. The 
original invariant encodes combinatorial information, and the multivariable gen- 
eralization will encode not only the same information but also more refined in- 
formation. For example, the chromatic polynomial of any tree with n vertices 
has chromatic polynomial x{x — 1)"^^. But not all trees have the same Xg{x). 
For example, if Ki^^ is the 4-star graph and is the path of order 4, Xk^ 3(2^) 
has a term XiXj for all i ^ j, but such a term is not present in Xp4.{x) . In fact, 
it is still an open question if X distinguishes trees, that is if Xti{x) ^ Xx2{x), 
whenever Ti and T2 are not isomorphic trees. 

A similar multivariable extension of the bad coloring polynomial is also natu- 
ral, especially given the importance of the latter because of its being equivalent 
to the Tutte polynomial. The following generalization of the bad coloring 
polynomial is also due to Stanley [Sta98] . 

Definition 13 Let G = iV,E) be a graph, let (j) : V ^ V = {1,2, . . .}, and let 
b{(j)) be the set of monochromatic edges in the coloring given by (j). Then the 
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symmetric function generalization of the bad coloring polynomial over indeter- 
minates xi,X2, ■ ■ ■ and t is 

XGix,t)= ^ + 

where the sum is over all possible colorings <j) of the graph G. 

Again, by setting Xi = 1 for 1 < i < A and xj = for j > A wc get the 
bad-coloring polynomial, and hence the Tutte polynomial. Therefore, Xg{x, t) 
is a multivariable generalization of the Tutte polynomial. 

There is another multivariable generalization of the Tutte polynomial that 
was developed independently and for very different reasons. This generalization 
is called the U-polynomial and is due to Noble and Welsh in [NW99]. 

Definition 14 Let G = {V, E) be a graph. Then the U-polynomial of G is 
C/G(a;,y)= ^a;„,---x„,(2/-l)l^l-'-(^), 

A<ZE 

where rii, . . . , are the numbers of vertices in the k different components of G 
restricted to A. 

Clearly, this is a generalization of the Tutte polynomial, as by setting Xi = 
{x — 1) for all i in Uaix^y) we get {x - 1)''("^)Tg(x, y). Note that the factor 

■ ■ ■ in every term keeps track of the number of vertices in the different 
components in A. This is a refinement of the rank-nullity generating-function 
definition of the Tutte polynomial where the factors a;'"("^)~'"(^) = 2;k(^)-«(G) 
each term keep track of the number of components in A. 

That Ug captures more combinatorial information from G than the Tutte 
polynomial can be seen by noting that Ug contains the matching generating 
polynomial, and thus the matching polynomial, as a specialization as well. 

Theorem 14 For any graph G, 

5(G;a;) = f/G(l,t,0,...,0,...,y=l). 

The U-polynomial has a deletion/contraction reduction relationship not in 

the class of graphs but in the class of weighted graphs. To see this, we turn to 
the W-polynomial also due to Noble and Welsh in [NW99]. A weighted graph 
consists of a graph G = {V, E), together with a weight function w : F — > Z+. 

If e is an edge of (G, uj) then (G \ e, uj) is the weighted graph obtained from 
(G, by deleting e and leaving ui unchanged. If e is not a loop, (G/e, w/e) is 
the weighted graph obtained from (G, w) by contracting e, that is deleting e and 
identifying its endpoints v, v' into a single vertex v" . The weight function oj/e 
is defined as ui/e{u) = U!{u) for ail u £V \ {v, v'} and Lo/e{v") = lv{v) + lo{v'). 
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Definition 15 Let (G, lu) be a weighted graph. The W-polynomial may be given 
recursively by the following rules. If e is an ordinary edge or a bridge, then 

Ife is a loop, then W(^q^i^){x, y) = yW(G\e,u){x, y). Otherwise, (G, lo) is En, the 
edgeless graph on n > vertices, with weights ai, . . . ,a„ and W(^E„.ui){x,y) = 

^ai * * ■ ^aji • 

That the resulting multivariate polynomial W is independent of the order in 
which the edges are deleted and contracted is proved in [NW99] . This can easily 
be done by induction on the number of edges once it is proved that the order in 
which you contract or delete edges in (G, lS) does not affect the weighted graph 
which you obtain. 

The U-polynomial is obtained from the W-polynomial by setting all weights 
equal to 1 and a proof that this definition is equivalent to Definition [14] can be 
found in [NW99] . Actually in [NW99] it is proved that W has a representation 
of the form 

W(G,..){x,y) = x.^---^c,{y-l)\^^-'-^^\ 

ACE 

where Ci, 1 < i < fc, is the total weight of the ith component of the weighted 
subgraph (A, w). 

Noble and Welsh [NW99] show that the symmetric function generalization 
of the bad coloring polynomial and the U-polynomial are equivalent in the fol- 
lowing sense. 

Theorem 15 For any graph G, the polynomials Ug and Xg determine each 
other in that if Po — ^ and pr — x"^ , then 

XG{x,t)^t\''\UG{x,^^,y = t + l). 

There is yet another polynomial, the polychromate, introduced originally by 
Brylawski in [BrySl], that is as general as Ug or Xq. Given a graph G and 
a partition vr of its vertices into non-empty blocks, let e(7r) be the number of 
edges with both ends in the same block of the partition. If t{tt) ~ (tii, . . . , ti^) 
is the type of partition tt, we denote by a;T-(^) the monomial Y[i=i ■ 

Definition 16 Let G be a graph. Then the polychromate XG{x,y) is 

XG{x,y) = ^y'^^'^^a^rW, 

where the sum is over all partitions ofV{G). 

We have the following theorem due to Sarmiento in [SarOO] but see [MN] for 
a different proof. 

Theorem 16 The polynomials UG{x,y) and XG{x,y) are equivalent. 
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The story doesn't end here. All three polynomials Ucix.y), XG{x,t) and 
XG{x,y) have natural extensions. For example the extension of the XG{x,t) 
replaces the t variable by countably infinitely many variables ti,t2,---, enu- 
merating not just the total number of monochromatic edges but the number of 
monochromatic edges of each color. It is defined as follows. 



where the sum is over all colorings ^ of G and is the set of monochromatic 
edges for which both end points have color i. By setting ti = t for alH > 1 we 
regain XG{x,t). 

For the other extension the reader is referred to [MN]. There it is also proved 
that all these extensions are equivalent. 

4.2 The Pcirametrized Tutte Polynomial 

The basic idea of a parametrized Tutte polynomial is to allow each edge of a 
graph to have four parameters (four ring values specific to that edge), which 
apply as the Tutte polynomial is computed via a deletion/contraction recursion. 
Which parameter is applied in a linear recursion reduction depends on whether 
the edge is deleted or contracted as an ordinary edge, or whether it is contracted 
as an isthmus or deleted as a loop. The difficulty lies in assuring that a well 
defined function, that is, one independent of the order of deletion/contraction, 
results. This requires a set of relations, coming from three very small graphs, 
to be satisfied. Interestingly, additional constraints are necessary for there to 
be a corank-nuUity expansion or even for the function to be multiplicative or a 
graph invariant, that is, equal on isomorphic graphs. 

The motivation for allowing edge-specific values for the deletion/contraction 
recursion comes from a number of applications where it is natural. This includes 
graphs with signed edges coming from knot theory, graphs with edge-specific fail- 
ure probability in network reliability, and graphs whose edges represent various 
interaction energies within a molecular lattice in statistical mechanics. While 
there is compelling motivation for allowing various edge parameters, the tech- 
nical details of a general theory are challenging. The two major works in this 
area are Zaslavsky [Zas92] and BoUobas and Riordan [BR99]. However, these 
two works take different approaches, which were subsequently reconciled with 
a mild generalization in [E-MT06], and for this reason we adopt the formal- 
ism of [E-MT06]. Bollobas and Riordan [BR99] also give a succinct historical 
overview of the development of these multivariable extensions. 

For the purposes of the following, we consider a class of graphs minor-closed 
if it is closed under the deletion of loops, the contraction of bridges, and the con- 
traction and deletion of ordinary edges; however we do not require closure under 
the deletion of bridges. Some formalism is necessary to handle the parameters. 

Definition 17 Let U be a class, and let R be a commutative ring. Then an R- 
parametrization of U consists of four parameter functions x,y,X,Y : U ^ R, 
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denoted e Xe,ye, Xg, Kg. 

Definition 18 Let U be an R-parametrized class, and let T he a minor-closed 
class of graphs with E{G) C U for all G € T. Then a parametrized Tutte 
polynomial on T is a function T : T —i- R which satisfies the following: T{G) = 
XeT(G/e) for any bridge e of G €T, and T{G) = YeT(G — e) for any loop e of 
G €T, and T(G) — y^T^G — e) + x^T^G/e) for any ordinary edge e. 

The following theorem gives the central result. Identity in Item [T] comes 
from requiring to be equal the two ways of carrying out deletion/contraction 
reductions on a graph on two vertices with two parallel edges ei and 62 having 
parameters {xei,yei, Xe^jYe-}. Similarly, identities in Items [2] and [s] come from 
considering the 0-graph and K3. Here again En is the edgeless graph on n 
vertices. 

Theorem 17 (The generalized Zaslavsky-Bollobas-Riordan theorem) 

Let R be a commutative ring, let T be a minor-closed class of graphs whose edge- 
sets are contained in an R-parametrized class U , and let ai, 02, . . . € R. Then 
there is a parametrized Tutte polynomial T onT with T(En) — an for all n with 
En G r z/ and only if the following identities are satisfied: 

1. Whenever ei and 62 appear together in a circuit of a k-component graph 
G e r, then a^ixe^Ye^ + yeiXe.^) = akixe^Ye^ + ye^XeJ. 

2. Whenever ei, ei and 63 appear together in a circuit of a k-component 

graph G G T, then akX^^{x^^Ye^ + ye^Xe^) = Ofc-'^ea (i^eiS^es + a^eiJ/ea)- 

3. Whenever Ci, 62 a,nd 63 are parallel to each other in a k-component graph 
G GT, then akYe^{XeiYe^ + ye^Xe^) — akYesO^e 

A most general parametrized Tutte polynomial, what possibly could be 
called the parametrized Tutte polynomial, might begin with the polynomial 
ring on independent variables {xe-,ye,Xe,Ye: e G C/} U {a^ : i ^ 1}. However, 
the resulting function is not technically a polynomial, in that it must take its 
values not in the polynomial ring, but has as R the polynomial ring modulo the 
ideal generated by the identities in Theorem [TTj 

The question also arises as to whether "the most general" parametrized Tutte 
polynomial should be multiplicative on disjoint unions and one point joint of 
graphs, as this introduces additional relations among the a^'s. This is because 
a parametrized Tutte polynomial is not necessarily multiplicative. A sufficient 
condition is the following. 

Proposition 4 Suppose T is a parametrized Tutte polynomial on a minor- 
closed class of graphs that contains at least one graph with k components for 
every k and that is closed under one-point unions and the removal of isolated 
vertices. Then T is multiplicative with respect to both disjoint unions and one- 
point joins if and only if the ai — T{Ei) is idempotent, and au = ai for all 
fc > 1. 
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BoUobas and Riordan [BR99] emphasize graph invariants, and hence require 
that the parametrization be a coloring of the graph. That is, graphs are edge- 
colored (not necessarily properly) , with edges of the same color having the same 
parameter sets. This enables consideration of parametrized Tutte polynomials 
that are invariants of colored graphs, but requires the following additional con- 
straints. For every ei € U, there are 62, 63 G C/ with ei 7^ 62 7^ 63 7^ ei such 

that Xq-^ — — "^es 5 2/ei — Ve^ — ' — -^e2 — and — — ^£3' 

Proofs of the above results and further details may be found it [Zas92,BR99, 
E-MT06] . We note that any relation between this Tutte polynomial generaliza- 
tion with its edge parameters, and the W- and U-polynomials of Subsection |4.1| 
with their vertex weights, has not yet been studied. 

Interestingly, although the parametrized Tutte polynomial has an activities 
expansion analogous to that of the classical Tutte polynomial, it does not nec- 
essarily have an analog of the rank- nullity formulation. However, under modest 
assumptions involving non-zero parameters and some inverses, the parametrized 
Tutte polynomial may be expressed in a rank-nullity form. This is fortunate, 
because significant results for the zeros of the chromatic and Tutte polynomial 
have arisen from such a multivariable realization. Examples may be found in 
Sokal [SokOla], Royle and Sokal [RS04], and Choe, Oxley, Sokal, and Wag- 
ner [COSW04]. 

4.3 The Generalized Transition Polynomial 

A number of state model polynomials, for example the circuit partition poly- 
nomials, Penrose polynomial, the Kauffman bracket for knots and links, and 
the transition polynomials of Jaeger [Jae90], that are not specializations of the 
Tutte polynomial, are specializations of the multivariable generalized transition 
polynomial of [E-MS02] which we describe here. This multivariable extension 
is a Hopf algebra map, which leads to structural identities that then inform its 
various specializations. The medial graph construction that relates the circuit 
partition polynomial and the classical Tutte polynomial extends to similarly 
relate the generalized transition polynomial and the parametrized Tutte poly- 
nomial when it has a rank- nullity formulation. 

The graphs here are Eulerian, although not necessarily connected, with 
loops and multiple edges allowed. A vertex state, or transition, is a choice of 
local reconfiguration of a graph at a vertex by pairing the half edges incident 
with that vertex. A graph state, or transition system, S{G), is the result of 
choosing a vertex state at each vertex of degree greater than 2, and hence is a 
union of disjoint cycles. We will write St{G) for the set of graph states of G, 
and throughout we assume weights have values in R, a commutative ring with 
unit. 

A skein relation for graphs is a formal sum of weighted vertex states, together 
with an evaluation of the terminal forms (the graph states). See [E-M98, E- 
MS02] for a detailed discussion of these concepts, which are appropriated from 
knot theory, in their most general form, and Yetter [Yet90] for a general theory 
of invariants given by linear recursion relations. A skein type, (or state model. 
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or transition) polynomial is one which is computed by repeated appHcations of 
skein relations. See Jaeger [Jae90] for a comprehensive treatment of these in 
the case of 4-regular graphs. 

For brevity, we elide technical details such as free loops and isomorphism 
classes of graphs with weight systems which may be found in [E-MS02]. 

Definition 19 Pair, vertex, and state weights: 

1. A pair weight is an association of a value p(e^,e^) in a unitary ring R 
to a pair of half edges incident with a vertex v in G. A weight system, 
W{G), of an Eulerian graph G is an assignment of a pair weight to every 
possible pair of adjacent half edges of G. 

2. The vertex state weight of a vertex state is Y[p{^vj e^) where the product 
is over the pairs of half edges comprising the vertex state. 

3. The state weight of a graph state S of a graph G with weight system W 

is Lu(S) = J3a;(u, S), where lo{v, S) is the vertex state weight of the vertex 
state at v in the graph state S, and where the product is over all vertices 



When A is an Eulerian subgraph of an Eulerian graph G with weight system 

W{G), then A inherits its weight system T'K(A) from G in the obvious way, 
with each pair of adjacent edges in A having the same pair weight as it has in 
G. When A is a graph resulting from locally replacing the vertex v by one of 
its nr^To'^ (2" ~ (2'i + 1)) vertex states, then all the pair weights are the same 
as they are in G, except that all the pairs of half edges adjacent to the newly 
formed vertices of degree 2 in A have pair weight equal to 1, the identity in R. 

The generalized transition polynomial N (G; W, x) has several formulation, 
and we give two of them, a linear recursion formula and a state model formula, 
here. 

Definition 20 The generalized transition polynomial, N{G;W,x), is defined 
recursively by repeatedly applying the skein relation 



at any vertex v of degree greater than 2. Here the Gi 's are the graphs that result 
from locally replacing a vertex v of degree 2n in G by one of its vertex states. The 
Pi 's are the vertex state weights. Repeated application of this relation reduces 
G to a weighted (formal) sum of disjoint unions of cycles, (the graph staies). 
These terminal forms are evaluated by identifying each cycle with the variable 
X, weighted by the product of the pair weights over all pairs of half edges in the 
cycle. 

Definition 21 The state model definition of the generalized transition polyno- 
mial is: 



ofG. 



NiG;W,x) = Y,PiN{Gi;W{Gi),x) 
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Note that vertex states commute, that is, if Guv results from choosing a 



vertex state at u, and then at v, we have Guv — Gvu- Thus, Definition 20 gives 
a well-defined function, and Definitions [20| and [21] are equivalent. 

Several of the polynomials we have already seen are specializations of this 
generalized transition polynomial. For example, if all the pair weights are 1, 
then the circuit partition polynomial for an unoriented Eulerian graph results. 
If G is an Eulerian digraph, and G is the underlying undirected graph with pair 
weights of 1 for pairs half edges corresponding to one inward and one outward 
oriented half edge of G and otherwise, then the oriented version of the circuit 
partition polynomial results. 

In the special case that G is 4-regular, the polynomial N (G; W, x) is es- 
sentially the same as the transition polynomial Q{G,A,t) of Jaeger [Jae90], 
where G is a 4-regular graph and A is a system of vertex state weights (rather 
than pair weights). If the vertex state weight in {G,A) is w, then define W{G) 
by letting the pair weights for each of the two pairs of edges determined by 
the state be ^/w. The two polynomials then just differ by a factor of a;, so 
N (G; W, x) = xQ{G, A, x), and here we retain vertices of degree 2 in the recur- 
sion while they are elided in [Jae90]. Thus N {G;W, x) gives a generalization of 
Jaeger's transition polynomials to all Eulerian graphs. 

Because Q{G,A,t) assimilates the original Martin polynomial for 4-regular 
graphs and digraphs, the Penrose polynomial and the Kauffman bracket of 
knot theory (see [Jae90]), and N {G;W,x) assimilates Q{G,A,t), we have 
that the Penrose polynomial and Kauffman bracket are also specializations of 
(G; W, x) . Specifically, if G is a planar graph with face 2-colored medial graph 
Gm, and we give a weight system to Gm by assigning a value of 1 to pairs of 
edges that either cross at a vertex or bound the same black face and otherwise, 
then N {Gm',W,x) = P{G;x). Similarly, if L is a link, and Gl is the signed, 
face 2-colored universe of L, then a weight system can be assigned to Gl so that 
N {Gl; W, + a'^) = (a^ + a-^)K[L] where K[L] is the Kauffman bracket of 
the link. 

Because of these specializations, the algebraic properties of the generalized 
transition polynomial are available to inform these other polynomials as well. In 
particular, N {G; W,x) is a Hopf algebra map from the freely generated (com- 
mutative) hereditary Hopf algebra of Eulerian graphs with weight systems to 
the binomial bialgebra R[x] (details may be found in [E-MS02]). This leads to 
two structural identities, the first from the comultiplication in the Hopf algebra, 
the second from the antipode. 

Theorem 18 Let G be an Eulerian graph. Then 

TV (G; ly, a; + y) - ^ TV (Ai; (Ai) ,x)N{A2;W {A2) , y) 

where the sum is over all ordered partitions ofG into two edge-disjoint Eulerian 
subgraphs Ai and A2, and 

N{G;W,-x)=N{C{G; W),x), 
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where C is the antipode C{G,W) = ^ (— l)'^' (Ai . . . A|p|), with the sum over 
all ordered partition P of G into \P\ edge-disjoint Eulerian subgraphs each with 
inherited weight system. Here N{G; W, x) is extended linearly over such formal 
sums. 

This type of Hopf algebraic structure has already been used to considerably 
extend the known combinatorial interpretations for evaluations of the Martin, 
Penrose, and Tutte polynomials implicitly by BoUobas [Bol02], and explicitly by 
Ellis-Monaghan and Sarmiento [E-M98,E-MS01, SarOl, E-M04a, E-M04b]. The 
first identity has been used to find combinatorial interpretations for the Martin 
polynomials for all integers, where this was previously only known for 2, -1, 
0, 1 in the oriented case, and 2, 0, 2 in the unoriented case. This then led to 
combinatorial interpretations for the Tutte polynomial (and its derivatives) of 
a planar graph for all integers along the line x = y, where previously 1, 3 were 
the only known non-trivial values. These results for the Tutte polynomial were 
mentioned in the preceding chapter and for the circuit partition polynomial in 



Subsection 3.4 The second identity has been used to determine combinatorial 
interpretations for the Penrose polynomial for all negative integers, where this 
was previously only known for positive integers. 



5 Two Applications 

Graph polynomials have a wide range of applications throughout many fields. 
We have already seen some examples of this with various applications of the 
classical Tutte polynomial in the previous chapter. Here we present two rep- 
resentative important applications (out of many possible) and show how they 
may be modeled by graph polynomials. 

5.1 DNA Sequencing 

We begin with string reconstruction, a problem that may be modeled by the 
interlace and circuit partition polynomials (and hence indirectly in special 
cases by the Tutte polynomial) . String reconstruction is the process of reassem- 
bling a long string of symbols from a set of its subsequences together with some 
(possibly incomplete, redundant, or corrupt) sequencing information. While 
we focus on DNA sequencing, which was original the motivation for the de- 
velopment of the interlace polynomial, the methods here apply to any string 
reconstruction problem. For example, fragmenting and reassembling messages 
is a common network protocol, and reconstruction techniques might be applied 
when the network protocol has been disrupted, yet the original message must 
be reassembled from the fragments. 

DNA sequences are typically too long to read at once with current labora- 
tory techniques, so researchers probe for shorter fragments (reads) of the strand. 
They then are faced with the difficulty of recovering the original long sequence 
from the resulting set of subsequences. DNA sequencing by hybridization is a 
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method of reconstructing the nucleotide sequence from a set of short substrings 
(see Waterman [Wat95] for an overview). The problem of determining the num- 
ber of possible reconstructions may be modeled using Eulerian digraphs, with 
a correct sequencing of the original strand corresponding to exactly one of the 
possible Eulerian circuits in the graph. The probability of correctly sequencing 
the original strand is thus the reciprocal of the total number of Euler circuits 
in the graph. 

The most basic (two-way repeats only) combinatorial model for DNA se- 
quencing by hybridization uses an Eulerian digraph with two incoming and 
two outgoing edges at each vertex (see Pevzner [Pev89] and Arratia, Martin, 
Reinert, Waterman [AMKW96]). The raw data consists of all subsequences of 
the DNA strand of a fixed length L, called the L-spectrum of the sequence. 
As L increases, the statistical probability that the beginning and end of the 
DNA strand are the same approaches zero, as does the likelihood of three or 
more repeats of the same pattern of length L or more in the strand (see Dyer, 
Frieze, and Suen [DFS94]). Thus, this model assumes that the only consid- 
eration in reconstructing the original sequence is the appearance of interlaced 
two-way repeats, that is, alternating patterns of length L or greater, for example, 
. . . ACTG . . . CTCT . . . ACTG . . . CTCT .... 

From the multiset (duplicates are allowed) of subsequences of length L, create 
a single vertex of the de Bruijn graph for each subsequence of length L-1 that 
appears in one of the subsequences. For example, if L = 4 and ACTG appears 
as a subsequence, create two vertices, one labeled ACT and one labeled CTG. 
Edges are directed from head to tail of a subsequence, e.g. there would be a 
directed edge from ACT to CTG labeled ACTG. If there is another subsequence 
ACTT, we do not create another vertex ACT, but rather draw an edge labeled 
ACTT from the vertex ACT to a new vertex labeled CTT. If, in the multiset 
of subsequences, ACTG appears twice, then we draw two edges from ACT to 
CTG. 

The beginning and end of the strand are identified to be represented by 

the same vertex, and, since by assumption no subsequence appears more than 
twice, the result is an Eulerian digraph of maximum degree 4. Tracing the 
original DNA sequence in this graph corresponds to an Eulerian circuit that 
starts at the vertex representing the beginning and end of the strand. All 
other possible sequences that could be (mis)reconstructed from the multiset of 
subsequences correspond to other Eulerian circuits in this graph. Thus (up to 
minor reductions for long repeats and forced subsequences) , finding the number 
of DNA sequences possible from a given multiset of subsequences corresponds 
to enumerating the Eulerian circuits in this directed graph. 

The generalized transition polynomial models this problem directly: when 
the pair weights are identically 1, it reduces to the circuit partition polynomial. 
This is a generating function for families of circuits in a graph, so the coefficient 
of X is the number of Eulerian circuits. The interlace polynomial informs the 
problem as follows. Consider an Eulerian circuit through the de Bruijn graph, 
which gives a sequence of the vertices visited in order. Now construct the inter- 
lace graph by placing a vertex for each symbol and an edge between symbols that 
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are interlaced (occur in alternation) in the sequence. The interlace polynomial 
of the interlace graph is then a translation of the circuit partition polynomial 
of the original de Bruijn graph, as in Theorem |12| where again the coefficient of 
X is the number of Eulerian circuits (see Arratia, BoUobas, Coppersmith, and 
Sorkin [ABCSOO], and Arratia, BoUobas, and Coppersmith [ABSOO, ABS04a]). 

One of the original motivating goals of Arratia, BoUobas, Coppersmith, and 
Sorkin [ABCSOO] was classifying Eulerian digraphs with a given number of Eule- 
rian circuits. The BEST theorem, a formula for the number of the circuits of an 
Eulerian graph in terms of its Kirchhoff matrix (see Fleischner [Fle91] for good 
exposition) gives only a tautological classification: the Eulerian digraphs with 
TO Eulerian circuits are those where BEST theorem formula gives to circuits. 
Critically, all of the above graph polynomials encode much more information 
than is available from the BEST theorem, and all of them are embedded in 
broader algebraic structures that provide tools for extracting information from 
them. Thus, they better serve the goal of seeking structural characterizations 
of graph classes with specified Eulerian circuit properties. 



5.2 The Potts Model of Statistical Mechanics 

Here we have an important physics model that remarkably was found to be 
exactly equivalent to the Tutte polynomial. 

Complex systems are networks in which very simple interactions at the mi- 
croscale level determine the macroscale behavior of the system. The Potts model 
of statistical mechanics models complex systems whose behaviors depend on 
nearest neighbor energy interactions. This model plays an important role in 
the theory of phase transitions and critical phenomena, and has applications as 
widely varied as magnetism, adsorption of gases on substrates, foam behaviors, 
and social demographics, with important biological examples including disease 
transmission, cell migration, tissue engulfment, diffusion across a membrane, 
and cell sorting. 

Central to the Potts model is the Hamiltonian, 

h{u}) ^ - J ^ 5{(Ji,(Jj) , 

{i,3}eE[G) 

a measure of the energy of the system. Here a spin, dj, at a vertex i, is a choice 
of condition (for example, healthy, infected or necrotic for a cell represented 
by the vertex). J is a measure of the interaction energy between neighboring 
vertices, w is a state of a graph G (that is, a fixed choice of spin at each vertex), 
and 5 is the Kronecker delta function. 

The Potts model partition function is the normalization factor for the Boltz- 
mann probability distribution. Systems such as the Potts model, following 
Boltzmann distribution laws, will have the number of states with a given en- 
ergy (Hamiltonian value) exponentially distributed. Thus, the probability of 
the system being in a particular state uj at temperature t is: 

exp (— {lo)) 



Pr(c^,/3) 



X;exp(-/3/i (tu))' 
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Here, the sum is over all possible states w of G, and 13=-^, where n = 
1.38 X 10~^^ joules/Kelvin is the Boltzmann constant. The parameter t is an 
important variable in the model, although it may not represent physical temper- 
ature, but some other measure of volatility relevant to the particular application 
(for example ease of disease transmission/reinfection). The denominator of this 
expression, P{G;q,P) = ^ exp (— (ro)), called the Potts model partition 
function, is the most critical, and difficult, part of the model. 

Remarkably, the Potts model partition function is equivalent to the Tutte 
polynomial: 

PiG;q,(3) - gMG)„k(G)|-/c(G)j. (^G; ^,«+ , 

where q is the number of possible spins, and v = exp(J/3) — 1. See Fortuin 
and Kasteleyn [FK72] for the nascent stages of this theory, later exposition in 
Tutte [Tut84], Biggs [Big96], BoUobas [Bol98], Welsh [Wel93], and surveys 
by Welsh and Merino [WMOO], and Beaudin, EUis-Monaghan, Pangborn and 
Shrock [BE-MPS]. 

One common extension of the Potts model involves allowing interaction en- 
ergies to depend on individual edges. With this, the Hamiltonian becomes 
h{uj) = J2eeE{G) JeS{cri, aj), where Jg is the interaction energy on the edge e. 
The partition function is then 

P{G)= Y: 'I'^^'Il-e, 

ACE(G) eeA 

where Ve = exp(/3Je) — 1 Again see Fortuin and Kasteleyn [FK72], and more 
recently Sokal [SokOO, SokOlb]. As we have seen in Subsection |4.2[ the Tutte 
polynomial has also been extended to parametrized Tutte functions that incor- 
porate edge weights. The generalized partition function given above satisfies the 
relations of Theorem |17| however, and thus is a special case of a parametrized 
Tutte function. 

This relationship between the Potts model partition function and the Tutte 
polynomial has led to a remarkable synergy between the fields, particularly for 
example in the areas of computational complexity and the zeros of the Tutte 
and chromatic polynomials. For overviews, see Welsh and Merino [WMOO], and 
Beaudin, EUis-Monaghan, Pangborn and Shrock [BE-MPS]. 



6 Conclusion 

There are a great many other graph polynomials equally interesting to those 
surveyed here, including for example the F-polynomials of Farrell, the Hosaya 
or Wiener polynomial, the clique/ independence and adjoint polynomials, etc. 
In particular, Farrell [Far79b] has a circuit cover polynomial (different from 



the circuit partition polynomial of Subsection 3.4 1 with noteworthy interrela- 



tions with the characteristic polynomial. Also, Chung and Graham developed 
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a 'Tutte-like' polynomial for directed graphs in [CG95]. The resultant cover 
polynomial is extended to a symmetric function generalization, like those in 
section [il] by Chow [Cho96]. Similarly, Courcelle [Cou08] and Traldi [Tra] 
have also very recently developed multivariable extensions of the interlace poly- 
nomial. Some surveys of graph polynomials with complementary coverage to 
this one include: Pathasarthy [Par89], Jaeger [Jae90], Farrell [Far93], Fiol 
[Fio97], Godsil [God84], Aigner [AigOl], Noy [NoyOB] and Levit and Man- 
drescu [LM05]. 
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